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Abstract 

We present a model of cosmic ray (CR) injection into the Galactic space based 
on recent 7-ray observations of supernova remnants (SNRs) and pulsar wind 
nebulae (PWNe) by the Fermi Large Area Telescope [Fermi) and imaging 
atmospheric Cherenkov telescopes (IACTs). Steady-state injection of nuclear 
particles and electrons (e~) from the Galactic ensemble of SNRs, and elec- 
trons and positrons (e + ) from the Galactic ensemble of PWNe are assumed, 
with their injection spectra inferred under guidance of 7-ray observations 
and recent development of evolution and emission models. The ensembles of 
SNRs and PWNe are assumed to share the same spatial distributions. As- 
sessment of possible secondary CR contribution from dense molecular clouds 
interacting with SNRs is also given. CR propagation in the interstellar space 
is handled by GALPROP. Different underlying source distribution models 
and Galaxy halo sizes are employed to estimate the systematic uncertainty 
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of the model. We show that this observation-based model reproduces the 
positron fraction e + / (e~ + e + ) and antiproton-to-proton ratio (p/p) reported 
by PAMELA and other previous missions reasonably well, without calling 
for any speculative sources. A discrepancy remains, however, between the 
total e~ + e + spectrum measured by Fermi and our model below ~ 20 GeV, 
for which the potential causes are discussed. Important quantities for Galac- 
tic CRs including their energy injection, average lifetime in the Galaxy, and 
mean gas density along their typical propagation path are also estimated. 

Keywords: ISM: cosmic rays, ISM: supernova remnants, ISM: clouds, 
7-rays: observations 



1. Introduction 

Recent observation of the positron fraction, e + /(e~ + e + ), by PAMELA 
[l[ shows an excess (referred to as the e + excess, or the excess in the positron 
fraction) relative to the prediction of a cosmic-ray (CR) propagation model 
jij in the energy range between 10 and 100 GeV. The model referenced in 
most analyses of the e + excess is GALPROP version which assumes e + 
to be produced along the propagation path by the Galactic CR protonffl 
having a power-law spectrum with a locally observed index of 2.75. A con- 
straint applied to these analyses is that the spectrum predicted for e~ + e + 



1 We refer to specific versions of GALPROP by the year of publication if not labeled in 
the literature. A detailed description of different versions of GALPROP can be found at 
http://GALPROP.stanford.edu, Various results obtained with GALPROP are reviewed 
in Strong et al. [3]. 

2 Contribution of alpha particles and heavy ions to pion production is included in 
"protons" by scaling the cross-section by an effective 'nuclear enhancement factor' of 
f.68 [4j. We note that Mori [H| derived a larger factor (1.845) for CR energy >^ 
5 GeV/nucleon and that the recent CR measurements report spectral hardening at ener- 
gies >~ 200 GeV/nucleon 0,0. The present analysis is primarily concerned with CRs 
below ~ 300 GeV and statistically dominated by gamma-rays below ~ 1 GeV. The nuclear 
enhancement factor relevant to the present analysis should be obtained in the above CR 
and gamma-ray energy ranges and is probably smaller than 1.845 (see @,@|). At present 
we assume it to be 1.68. 

3 The known gamma-ray producing particle processes which do not go through neu- 
tral pions (e.g., if — > 77 and direct photon production processes) are also included but 
contribute less than 1 % in the present energy range. 
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or lepton^] agrees with that measured by Fermi jsl, [h| and the e~ spectrum 
measured by PAMELA [11]. Various additional sources of e + have been pro- 
posed to account for the e + excess, including pulsars (PSRs) and pulsar wind 



nebulae (PWNe) [e.g. Il2l. |13. Il4. Il5l. Il6l . Il7l . and references therein]; super- 
nova remnants (SNRs) [e.g. 



la . Il9l |20L and references therein] ; propagation 



effects [e.g., |2l|, |22j, and references therein]; and dark matter (DM) annihi- 



16l . |24J . and references therein] . In some references, 



lation or decay [e.g., C 
the e + excess is discussed together with a bump in the CR lepton spectrum 
claimed by ATIC (25|. T his bump has not been observed by Abdo et al. 



Ackermann et al. 10 in the latest Fermi measurements. Hence we will 



not consider the ATIC bump in this paper. 

Measurements of the antiproton-to- proton ratio (p/p) have recently been 
extended to ~ 100 GeV by PAMELA [26|, |27j. The reference model used in 



the analysis of the ratio is GALPROP98b which predicts p to be produced^] in 
the same inter-stellar matter (ISM) by the same Galactic CRs as for e + [e.g., 



30j,|3l|, and references therein]. Some early measurements of the ratio at lower 
energies at ~ 10 GeV gave higher values than the GALPROP98b prediction, 
and possible contribution from the DM annihilation have been discussed by 
Bergstrom et al. 32[ and in references given in Moskalenko et al. 



30 . The 



new PAMELA measurement agrees well with a recent GALPROP version 



labeled as DC in Moskalenko et al. [30j. Hence we will not consider possible 



contribution to CR p from the DM annihilation in this paper. 

In SNRs, 7, e~, e + , and p are produced in the pp interactions as secondary 
particles which can, in principle, account for the PAMELA e + excess. How- 
ever, the yields of tc + , tt~ and n° are predicted to be approximately equal for 
proton kinetic energies greater than ~ 10 GeV, hence the spectrum of e + , 
a daughter of tc + decay, is tightly constrained by the observed gamma-ray 
spectrum of 7i~° origin. We survey below SNR-based scenarios proposed to 
account for the e + excess. In the work by Blasi [18|, e + are assumed to be re- 



accelerated to enhance the positron fraction. The lepton spectrum will then 



4 We refer to e~ and e + collectively as "leptons" later in this paper. We note "leptons" 
include nominally muons, tau particles and neutrinos. 

5 We note that anti-neutrons are predicted to be produced about equally to or substan- 
tially more than p in the high energy pp interaction dependent upon the iso-spin nature 
of diquark pairs produced when the QCD color string breaks [11] ■ They decay to p with 
the lifetime of n. In GALPROP, the p inclusive cross section by Tan & Ng [29| has been 
doubled to include them. 
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become harder as energy increases until radiative cooling takes over in the 
multi TeV energy band. All charged CRs including p will also show similar 
spectral hardening [33;]. Fujita et al. 19 1 have assumed that e + are produced 



in dense clouds by nuclear CRs accelerated when the Local Bubble exploded. 
The dense clouds existing at that time are assumed to have been destroyed 
by now. In this scenario, we expect to find anisotropy in the arrival direction 
of high energy charged CRs and in the pionic gamma-ray emissivity at local 
molecular clouds; this will be tested in future Fermi observations. Ahlers et 
al. 
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have calculated the e + spectrum at Earth based on proton spectra 
deduced from the presumed pionic gamma-ray spectra of TeV SNRs. These 
SNRs, however, do not represent the Galactic ensemble of SNRs as seen in 



recent Fermi observations [34J, |35|, |36|, |37|, |38j. Katz et al. [21] claims that 
the e + excess comes out naturally if the radiative loss time is comparable to 
the propagation time (~ 10 7 yr) for e + of energy ~ 30 GeV. The authors 
have assumed a non-standard propagation yet to be confirmed. Stawarz et 
al. [22J note that the rollover in the high energy lepton spectrum @, B So) 
can be explained by the Klein-Nishina effect in the radiative loss. They then 
note that the observed e + excess can be reproduced only if CRs propagate 
through high density regions (n > 80 cm~ 3 ) and secondary positrons pass 
through regions where star light density is extremely high (energy density 
> 300 eV cm" 3 ). 

Another group of papers claim that one or a few PWNe within a few 



100 pc of the solar system can account for the e + excess [e.g. |12|, |l3j, Il6l. Il4. 
and references cited therein]. The PWNe are assumed to have accumulated 
high energy leptons over a few x 10 kyr and have released them in an interval 
shorter than the propagation time such that the highest-energy particles are 
now reaching Earth. We note that local PWNe and SNRs have been consid- 
ered as p ossible sources of CR positrons prior to the PAMELA experiments 
[e.g., |4l|, |42|, |43|. and references therein]. According to the recent theoretical 



studies [e.g., |44j, energetic leptons can be injected impulsively into Galactic 
space only when some condition is met in a short epoch in their life, and the 
probability of such a rare impulsive injection taking place right now from one 
of very few nearby PWNe is low. On the observational side, HESS has found 
several PWNe whose ages exceed ~ 30 kyr as discussed by Aharonian et al. 
45j, Gallant 46], Funk 47| and in references given in Mattana et al. 48 



Fermi has found GeV emission from Vela-X 49[ . The leptons responsible for 



the GeV emission have energy of a few xlO GeV and reside in the halo of 



the PWN adjacent to the region where TeV emission has been found [49j, |50 
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These observations in the GeV-TeV band as well as recent theoretical studies 



44l . l5ll l52l |53| suggest CR leptons are released slowly from PWNe allowing 



a larger ensemble of Galactic PWNe than those within a few 100 pc of Earth 
to contribute to the CR lepton spectrum. 
Grasso et al. |] 



Malyshev et al. 14] considered a distribution of pulsars 



and PWNe within a few kpc of the solar system to be responsible for the 
PAMELA e + excess. Grasso et al. 16[ have fitted the positron fraction well 
by assuming an e~ spectrum softer than that used in GALPROPv44_500180, 



the reference widely used in this kind of analyses [31]. In this scenario, the e 



flux at higher e nerg ies are dominated by unidentified nearby PWNe. Recently 
Delahaye et al. [54| studied possible energy ranges that nearby PSRs, PWNe, 
and SNRs contributed to the CR electrons and positrons at Earth. Their 
contributions are predicted to give a flat positron fraction at energies greater 
than 10 GeV. 

has associated the e + excess to an impulsive injection of e + by a 



Ioka 55 



presumed historic GRB that created the Local Bubble. Possible association 



of the Local Bubble with a GRB has been suggested earlier [e.g., |56j, [57 
Perrot & Grenier [56} have predicted the CR e~ spectrum at Earth will be be 
harder at higher energies (E > 100 GeV) in such a scenario and emissivities 
of gamma-rays at local molecular clouds will show directional dependence at 
the ~ 50 % level. Future Fermi observations of gamma-ray emission from 
local clouds will detect such a large anisotropy if it exists. 

Numerous publications have attempted to associate the e + excess with 
annihilation and/or decay of dark matter particles [e.g. 



23, 16, 24, and ref- 



erences therein], in which a number of interesting possibilities are proposed. 
In this work, however, we will not discuss any dark matter related scenarios. 

All publications surveyed above attempted to reproduce the e + excess by 
interpreting a number of known local objects as CR sources using adjust- 
to-fit CR injection spectra and/or adopting non-standard CR propagation 
processes. Although we cannot rule out all such possibilities, it is important 
to study how recent gamma-ray observations constrain the positron fraction 
and the p/p ratio within the conventional framework. Another important 
issue with these publications is that spectra of 7, p, e~ , e + and p are analyzed 
more-or-less independently. If these CRs come from SNRs and PWNe, or are 
produced in nuclear interactions with molecular clouds and ISM gas, there 
will be strong correlations among their spectra. In particular, the CR 
spectrum should not be treated independently of the CR p, e + and p spectra. 
We put these data on the coherent platform of GALPROP by assuming the 
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Galactic CRs are in their steady states. 

This work focuses on the Galactic CR ratios (e + fraction and p/p ra- 
tio) and use published gamma-ray observation results as constraints and/or 
consistency check for our calculations. All Galactic CRs are assumed to be 
injected by the ensemble of SNRs and PWNe, or produced through their in- 
teractions with the interstellar gas in the Galaxy, or from clouds interacting 
with SNRs. The CR injection processes are all assumed to be in steady- 
state. The adopted injection spectrum of protons from SNRs is constrained 
by the available observations of CR proton flux at Earth, and is similar to 
those assumed in recent Fermi analyses of the diffuse Galactic gamma-ray 



emission 



58l l59l |60| . The injection spectrum of e + and e~ from PWNe is 
deduced from currently available gamma-ray observations of PWNe by Fermi 
and HESS, guided by recent theoretical spectral evolution models. The CR 
propagation process is calculated within the robust GALPROP framework. 
The calculated local CR fluxes are renormalized to those observed at Earth, 
which determines the CR injection luminosities. We purposely remain blind 
to the measured positron fraction as well as the measured p/p ratio until the 
calculated results are compared with the measurements. 

The major assumptions made in our model are summarized below: 

(i) CRs are injected from the Galactic ensembles of SNRs (primary parti- 
cles: p, e~; secondary particles e~, e + , p) and PWNe (primary particles: 
e~,e+). 

(ii) The primary p and e~ injected from the Galactic ensemble of SNRs have 
a common steady-state spectral shape with a fixed e/p ratio, except 
near the maximum energy where the e~ injection spectrum is corrected 
for radiative energy loss. 

(iii) The primary e~ and e + injected from the Galactic ensemble of PWNe 
have a common steady-state spectral shape and injection luminosity. 

(iv) The primary p spectrum injected from SNRs is constrained by the 
observed CR proton spectrum at Earth. 

(v) The secondary CR spectra injected from SNRs interacting with molec- 
ular clouds are calculated assuming the p injection spectrum obtained 
in (4) at the interacting sites. We use the pp — > e + ,e~ cross-sections 



modeled by Kamae et al. |6l|, |62| and the pp — >• p cross-section from 



Tan & Ng [29|, [63l |64 



(vi) The primary CR spectrum injected from PWNe is deduced from gamma- 
ray observations by Fermi and ACTs, and PWN evolution models, e.g. 
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by Tanaka & Takahara 52| and Zhang et al. 51 



(vii) SNRs and PWNe are distributed identically and continuously in the 
Galaxy using the parametric recipe implemented in GALPROP. Two 
different sets of parameters are used to test the robustness of the cal- 
culated results. 

(viii) The propagation, interaction and energy loss of CRs from the sources to 
Earth are calculated using GALPROP with three different Galaxy halo 
heights (2, 4 and 10 kpc). Propagation parameters in the GALPROP 
input, such as the diffusion coefficient and re- acceleration strength 
(Alfven velocity), are adjusted to reproduce the observed CR proton 
and lepton spectra, and the B/C ratio at Earth. 

This paper is organized as follows: We describe the three CR source 
classes assumed here in section |2j The two spatial distribution models of 
SNR/PWN in the Galaxy are described in section[3J CR injection from SNRs 
and PWNe are discussed in sections @] and 0, respectively. In section EJ the 
calculated positron fraction and the p/p ratio are compared with observations 
and discussion is given thereon. We also present quantitative calculations on 
energetics and live-times of Galactic CRs in the section. Finally, the results 
are concluded in section [7] with some comments on future prospects. 



2. Three Classes of Cosmic-ray Sources in the Galaxy 

The model consists of three classes of sources, each responsible for inject- 
ing a specific set of CR species into the interstellar space. We use GALPROP, 
whenever possible, to calculate CR injection and propagation. In some spe- 
cial cases, however, these are calculated using a simple propagation program 
specifically developed to accommodate, for example, non-power-law injection 
spectra and source distribution without cylindrical symmetry with respect 
to the Galactic Center. This will be discussed in more details in the relevant 
sections. 

The three classes of CR injection sources include: 

'SNR-propagation': This class includes the primary CRs (p and e~) in- 
jected from the Galactic ensemble of SNRs, the secondary CRs (e~, e + , 
p) and diffuse gamma-rays (7r° decay, bremsstrahlung and IC) produced 
along the propagation path of the primaries. 

'PWN-propagation': This includes the primary e~ and e + injected from 
the Galactic ensemble of PWNe, and the diffuse gamma-rays produced 
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along their propagation path via bremsstrahlung and inverse Compton 
scattering (IC). 

'SNR-cloud interaction': We assume that a fraction of CRs accelerated at 
SNRs (especially those in their middle-age) are interacting with dense 
local clouds (e.g., W51C, W44, IC 443, and W28). These interaction 
sites contribute through hadronic interactions to the injection of sec- 
ondary particles, including e +//_ , p, and gamma-rays. Propagation of 
e + /~ and p from this source class has been calculated with a simple 
propagation code described in Appendix. 

We use the conventional 2-dimensional propagation mode in GALPROP, 
for which the Galaxy boundary possesses cylindrical symmetry; we adopt a 
propagation volume bounded by R max = 30 kpc and Z max = 4 kpc (the halo 
height) for our 'Default' model. In addition, we also construct three alter- 
native models to check the robustness of our results to the assumptions made 
for the source distribution and the Galaxy halo size. The first alternative 
model adopts a different parameter set for the spatial distribution, while the 
second and third have the halo height modified to a smaller Z max = 2 and a 
larger Z max = 10 kpc respectively. We note that our results are checked to 
be insensitive to the change of the radial boundary R max by ±10 kpc. Prop- 
agation parameters such as the diffusion coefficient are adjusted for models 
with different halo heights in order to reproduce the observed B/C ratio. 

We also estimate, in the context of steady-state CR injection in the 
Milky Way, the significance of the injection of CR positrons from nearby 
sources (within a few hundred parsec from Earth) to the reproduction of 
the PAMELA positron fraction. Treatment of this CR contribution by these 
'local' sources in the steady-state regime will be discussed in section |6j 

3. Spatial Distribution of SNRs and PWNe 

We study two spatial distribution models of CR source, parameterized in 
the form of: 

Respectively, R and Z are the Galacto-centric radius and the distance from 
the Galactic Plane in kpc; R & = 8.5 kpc. The default model adopts a, 
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Figure 1 Probability density of source distribution plotted as a function of 
distance from Earth: the parametric model in GALPROP with the default 
(solid) and alternative model (dash) parameters. Corresponding curve de- 



rived from the 3D gas distribution obtained by Nakanishi & Sofue 66|, |67 
also shown dash-dotted line. 



is 



(3 = 1.25, 3.56 based on the pulsar distribution [31], while the alternative 
model fits to a SNR distribution model described in Case & Bhattacharya 



65( and uses a, (3 = 1.69, 3.33. The two source distribution models are 
shown in Fig. [1] as a function of distance from Earth. Also overlaid in the 
figure for reference purposes is a corresponding curve derived from the 3D 



gas distribution obtained by Nakanishi & Sofue [66j, [67 



4. Cosmic Ray Injection from Supernova Remnants 

4-1. Injection of Protons and Primary Electrons 

The standard theory of particle acceleration at SNRs through diffusive 
shock acceleration (DSA) predicts that the accelerated particles share a com- 
mon non-thermal spectral shape when gas density is not high (n < a few 
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10° 10 1 10 2 10 3 

E [GeV] 



Figure 2 The assumed steady-state injection spectrum of protons (solid) and 
primary electrons (dashed) for the 'SNR-propagation' source class. An ex- 
ponential cut-off at the maximum energy estimated for an age of 50 kyr is 
applied to the electron spectrum shown. 



cm 3 ), apart from a high-energy cutoff or rollover for electrons around their 



maximum energy limited by radiative loss [e.g., |68|, |69[. As a result, we 
adopt a common spectral shape for the SNR-injected protons and primary 
electrons, except near the maximum energy where radiative loss introduces 
a spectral cutoff for the e~ spectrum. We note that little is known about 
the cutoff at the highest energy theoretically and should be readjusted when 
observational data become available. 

We use a broken power-law injection spectrum for the SNR-injected pri- 
mary CRs. The spectral parameters and injection luminosity of protons are 
chosen such that the propagated spectrum fits well with the observed CR 
proton flux at Earth. The e/p ratio, K ep , is chosen such that the calculated 
total e~ + e + flux, which are mainly contributed by the 'SNR-propagation' 
and 'PWN-propagation' classes, reproduce the observed CR lepton spectrum 
measured by Fermi jsl, [To| in the 50 — 150 GeV energy band (more will be 
discussed in section [5]). The adopted steady-state proton and electron injec- 
tion spectrum are shown in Fig. [2] by a thick solid line and dash-dotted line 
respectively, and the relevant parameters are summarized in Table [TJ 
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Table 1. Summary of Model Parameters 



Parameter 



Value 



Description 



GALPROP Inputs 



B 

Z m . 
D 

VA 



x exp 



(M)xexp(-f) ,G 



50 

2, 4, 10 kpc 
(2.9 , 5.8, 10.0) x 10 28 /3 ( 
30 km s" 1 



_R_\C 

4GV / 



B-field model in ISM a 

Galactic halo height 
Diffusion coefficient b 
Alfven wave velocity in ISM 



Galactic Ensemble of Supernova Remnants 



Lp, >0.1GcV 


6.42 x 10 40 erg s" 1 


Total luminosity of proton 


L e , >0.1GcV 


1.23 x 10 39 erg s" 1 


Total luminosity of primary e - 


E-ep 


0.02 


The e/p ratio 


E b 


7 GeV 


Spectral break energy 


7i 


2.0 


Spectral index below E b 


72 


2.5 


Spectral index above E b 


n x W Pt >igov 


5 x 10 51 erg cm~ 3 


See footnote a in Table [2] 


Galactic Ensemble of Pulsar Wind Nebulae c d 


Le, >0.1GcV 


(2.38 - 5.72) x 10 37 erg s" 1 


Total luminosity of e + + e~ 


E b 


100 - 500 GeV 


Spectral break energy 


71 


1.0-1.7 


Spectral index below Eb 


72 


2.8 - 3.3 


Spectral index above Eb 



a R and Z are Galactocentric distances in kpc. 

b Parameter values for Z max = 2,4 and 10 kpc, respectively. 

c Quoted injection luminosities and K ep correspond to the default setup with Z n 
-Rmax = 30 kpc, and the default spatial source distribution. 



4 kpc, 



d Ranges of spectral parameters given represent the sampled possible intervals, which are as- 
sumed basing on gamma-ray observations of PWNe. The injection luminosity adjusts accordingly 
with the spectrum, with its out-coming range listed here. 
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According to Yamazaki et al. [70[ , the maximum electron energy is limited 
to ~ 14 TeV for SNRs older than 1 kyr, and ~ 7 TeV for those 10 kyr old; 
it then drops to ~ 0.5 TeV for ages of a few x 100 kyr. We calculate the 
loss-limited E max in accordance with Yamazaki et al. [70J (eqn. 2 and 10), 
assuming = 10 /iG, r = 4, h = E 5 i = v^g = n = 1, and a hydro solution 
for a typical SNR lifetime of 50 kyr within the Sedov phase, such that: 

B „„. 13 ,TeVx(J^-)"x(J^)" (2) 

Beyond 50 kyr, the electron injection luminosity is expected to be too low to 
influence the time-averaged injection; hence, it is neglected in this study. The 
evolution of -E max of the continuously injected e~ from the Galactic ensemble 
of SNRs results in a smooth rollover of the calculated local e~ spectrum from 
the 'SNR-propagation' class at around 1 TeV. To determine the shape of the 
rollover, we use the simple propagation code described in the Appendix. 
This rollover is then applied to the propagated primary e~ spectrum first 
calculated using GALPROP without the incorporation of E m3iX . More details 
of the calculation of the rollover shape is given in the Appendix. 

4-2. Secondary Particles Contribution from SNRs 

Secondary CRs are produced when the shock-accelerated nuclei interact 



with the surrounding gas. Blasi [18] proposed that secondary positrons pro- 
duced (and re-accelerated) in old SNRs may be responsible for the rise of the 
positron fraction observed by PAMELA. Meanwhile, antiprotons should also 



be produced and accelerated through an identical mechanism [33| . After the 
first year of science operation of the Fermi mission, high-energy gamma-ray 
emission have been detected with high statistical significance from 6 Galac- 
tic SNRs, and all of them have TeV counterparts observed by ground-based 
ACTs. Among them, 5 are middle-aged SNRs known to be interacting with 
local molecular clouds. The ir g amma-ray spectra are best explained by the 



decay of 7r° mesons [34j, |35|, |3g, |37|, |38| , with the underlying proton spectra 



at the interaction sites found to follow broken power-laws. A summary of 
the derived proton spectra for these gamma-ray SNRs is given in Table HJ 
The gamma-ray fluxes are found to be high, probably implying that the 
secondary CRs produced and injected by these interaction sites can be ap- 
preciable. However, the post-break indices of the deduced broken power-law 
proton spectra are generally soft (> 2.7), and particle acceleration in these in- 
teraction sites is probably inefficient beyond the break energies found around 
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~ 10 GeV (see Malkov et al. [71( for possible theoretical interpretation). We 
try to estimate in our model the contribution of secondary CRs from these 
interaction sites to the local fluxes (the 'SNR-cloud interaction' class), and 
deduce an upper-limit for the total number of these sites in our Galaxy using 
constraints from updated CR measurements. 

The injection luminosity of these secondary CRs are determined by the 
underlying proton spectrum and the product n x W p , where n and W p are 
the average gas density and the total kinetic energy of CR protons, within 
each interaction site. With only a relatively small sample of this source class 
available, our knowledge on these quantities is very limited. With reference 
to gamma-ray observations of SNRs interacting with clouds (Table [5]), we 
adopt a typical n x W p (T p > 1 GeV) of 5 x 10 51 erg cm -3 for each site. For 
this 'SNR-cloud interaction' source class, we use the same spatial distribu- 
tion (Section [3]), underlying proton spectrum, and propagation parameters 
as those adopted for the 'SNR-propagation' source class. Instead of using 
GALPROP, we use the simple propagation code described in the Appendix 
to propagate the secondary particles produced, such that the injection lumi- 
nosities per interaction site and the non-power-law injection spectra can be 
accommodated. As has been noted before, neither re-acceleration effect nor 
hadronic interaction along the CR propagation path are considered for these 
secondary CRs. 

5. Spectra of Electron and Positron Injected by PWNe 

Electrons and positrons are accelerated in the magnetosphere of a pul- 
sar, injected to its associated PWN, and then accelerated again, presumably 
through diffusive shock acceleration across a relativistic termination shock. 
Electromagnetic spectra in the keV and MeV bands are mostly believed to 
originate from synchrotron radiation by the accelerated leptons, while those 
in the GeV and TeV bands are dominated by IC scattering of the synchrotron 
as well as the background photon fields by the same lepton population. Un- 
til recently, the Crab nebula was the only PWN detected clearly in both the 
GeV and TeV bands, and little was known about the evolution of gamma-ray 
emission from PWNe. In the past few years, however, HESS has detected 



emission 



several PWN e [45l . l46l. l47j| . shedding light to the evolution of their broadband 



72, 48J. Mattana et al. 48] list 14 TeV PWNe of which Fermi has 



firmly detected GeV emission from Crab, Vela X and MSH 15-52 |73|, |49|, [74 



Of particular importance is that HESS has detected TeV emission from 3 
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Table 2. SNR-cloud interaction sites observed by Fermi and ACTs 



SNR 


Age 


Distance 


Gas dens x CR a 


PL or Broken- 


PL b 


Refs. c 


Name 


kyr 


kpc 


n x W p [erg cm~ 3 ] 


E br [GeV] 




OLHE 




W49B 


4.0 


7.5 


1.1 x 10 52 


4.0 


2.0 


2.7 


d 


W44 


20 


3.0 


6 x 10 51 


8.0 


1.74 


3.7 


e 


W51C 


30 


6.0 


5.2 x 10 51 


15 


1.5 


2.9 


f 


IC 443 


30 


1.5 


6.7 x 10 50 


69 


2.09 


2.87 


g 


W28 (North) 


40 


2.0 


1.3 x 10 51 


2.0 


1.70 


2.70 


h 



a Product of the average gas density and the total kinetic energy of CR protons in the SNR- 
cloud interaction site, assuming that the observed gamma-rays luminosities have a pion-decay 
dominated origin. Errors associated with these estimations are likely to be within a factor 
less than 10 above and below. 

b Spectrum of protons interacting with dense clouds at SNRs for hadronic scenarios. PL 
stands for a power-law spectrum. 

c References for ages and distances are found in the literatures cited below. 

d Abdo et al. dH 

e Abdo et al. [H 

f Abdo et al. 0, Fiasson et al. [75| 

sAbdo et al. 0, Albert et al. 0, Acciari et al. [z3 

h Abdo et al. [37|, Aharonian et al. [zi] 
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PWNe whose ages are estimated to be greater than 50 kyr. This suggests 
that PWNe retain TeV leptons up to or longer than ~ 10 5 yr. 

A number of recent studies have attempted quite successfully to construct 



dynamical evolution models of non-thermal emission of PWNe [5JJ, [15|, [52 , 



791 . |53| |. Some of these models are applied to explain mult i- wavelength spec- 
tra of a number of GeV or TeV PWNe (e.g. Crab Nebula, MSH 15-52, 
HESS J1825-137, Vela X, PWN GO.9+0.1 and G338.3+0.0) detected by 
ACTs and/or Fermi LAT, and obtain reasonably satisfactory results. The 
accumulated particle spectra derived from essentially all of these PWNe can 
be described by broken power-laws with the break energies lying around 
100 — 500 GeV. Tanaka & Takahara [52j used the broadband spectrum of 
the Crab Nebula to calibrate their evolution model, and found that a fast 
decaying magnetic field (~ t~ L5 ) and a slowly evolving spectrum of the ac- 
cumulated leptons (except for the first 1000 yr or so when B is still large) 
can possibly explain the trend of increasing TeV gamma-rays to X-ray ra- 
tio with age, as implied by recent multi-wavelength observations of PWNe. 
This weak time-dependence of the accumulated particle spectrum in PWNe 



has also been predicted by some other similar evolution model [e.g., [79j. In 
Table [3j we list out the spectral parameters for the accumulated particles 
predicted by a number of successful spectral evolution models for 5 repre- 
sentative gamma-ray PWNe observed in the TeV band. Among them, 4 of 
which also have firm detections by Fermi LAT. For all models, we choose the 
parameters for t > 1000 yr when the particle spectra become stable and are 
subjected to only weak time-dependence for the rest of the PWN lifetime. 

The actual process through which the accumulated energetic particles are 
released from the PWNe into the Galactic space is known very little, mainly 
because it is difficult to follow the evolution of a PWN lifetime through ob- 
servation, and also that the available samples of Galactic PWNe at different 
ages are only a handful. In the real situation, the real-time spectral shape 
of the 'escaped' particles may differ from that of the 'snapshot', effective 
spectra of the accumulated particles inside the PWNe directly derived by 
multi-wavelength modeling. However, the accelerated particles are expected 
to be completely released into the ISM within ~ 50 kyr, the typical lifetime 
of PWNe. As already mentioned above, spectral evolution models for PWNe 
have shown that time-evolution of the accumulated lepton spectrum inside 
the nebulae is expected to be very slow only for most part of a PWN life- 
time, due to a fast-decaying magnetic field. In the context of a straightly 
steady-state picture of our Galaxy, which we consider in this study, this subtle 
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Table 3. Accumulated lepton spectra in 5 PWNe detected by Fermi and 

ACTs 



Source 


Age 


Distance 


Firm detection 


Spectral parameters 


Refs. 


Name 


kyr 


kpc 


ACTs 


Fermi 


E br [GeV] 


OLLE 


OiHE 




Crab 


1.0 


2.0 


Y 


Y 


300 


1.5 


3.1 


b 


MSH 15-52 


1.7 


5.0 


Y 


Y 


460 


1.5 


2.9 


c 


HESS J1825-137 


21 


3.9 


Y 


N 


120 


1.0 


3.0 


d 


GO.9+0.1 


6.5 


8.5 


Y 


N 


80 


1.7 


3.3 


c 


G338.3-0.0 


4.5 


10.0 


Y 


Y 


260 


1.0 


2.8 


r 



a Effective spectra of the accelerated and accumulated electrons and positrons inside the 
PWNe, approximated by broken power-laws. For each source, if multiple references are 
given, the spectral model makes reference to the first listed. 



b Tanaka & Takahara [52j, Abdo et al. [73 1 
c Abdo et al. [71], Aharonian et al. [80j ] 

d Zhang et al. [51], Aharonian et al. 81, 82|, Grondin et al. [8£ 
e Qiao et al. [H] 

f Fang & Zhang 0, Slane et al. [H| 
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Figure 3 Spectra of accumulated lepton populations deduced from broadband 
emission modelings for the 5 PWNe listed in Table [3] (thin lines). All spectra 
are normalized to each other at 10 GeV for comparison. The yellow band 
shows the sampled possible range of steady-state injection spectral shape for 
leptons injected by the 'PWN-propagation' source class. The 'mean' spectral 
shape from the sampled band is shown by the thick line. 



time-dependence within such short time-scale is unimportant. We hence con- 
sider it a reasonable assumption that the snapshot momentum distribution of 
the accumulated particles can approximate to the first order the steady-state 
spectral shape of the e + /~ pairs being injected into the ISM by the Galactic 
PWN ensemble. We note that a similar prediction has been made before the 



PAMELA positron fraction measurement by Zhang & Cheng [86 



To determine this effective spectral shape, we refer to the measurement 
data and fitted parameters for the PWN sample given in Table[3]for guidance. 
We adopt a broken power-law distribution, and allow the spectral parameters 
to vary over conservative ranges spanned by the underlying particle spectra 
found from these 5 gamma-ray PWNe. Specifically, we vary the break energy 
over the 100 — 500 GeV band, and the pre-break and post-break power-law 
indices are sampled over values of 1.0 — 1.7 and 2.8 — 3.3, respectively. In 
Fig. [31 we show the underlying spectra derived for the 5 gamma-ray PWNe 
with thin lines. The range of steady-state injection spectrum we consider for 
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the PWN-propagation class is shown with the yellow shaded band. The thick 
line is the 'mean' spectral shape derived from this range of possible injection 
spectral shapes, which we will use as the reference spectrum (hereafter the 
'mean spectrum'). The steady-state CR injection luminosity for the Galactic 
PWN ensemble is determined jointly with the K ep parameter of the 'SNR- 
propagation' class by normalizing to the absolute lepton flux and reproduce 
the e + + e~ spectrum in the 50 — 150 GeV band measured by Fermi 0, Ho] ]. 
This procedure is sufficient to uniquely determine the absolute normalizations 
for both the 'SNR-propagation' and 'PWN-propagation' source classes. In 
the results of this study (Fig. [6j [7] and [8]) , we will use a gray band and a thick 
solid line to show the models which adopt the range of possible spectrum and 
the mean spectrum for the PWN-propagation CR injection class, respectively. 
The gray band will be considered as the systematic uncertainty of our result 
due to the uncertainty of the steady-state injection spectrum of the Galactic 
PWN ensemble. 

6. Results and Discussion 

The primary and secondary CR spectra obtained for the three source 
classes described in section [2] are summed and compared with the collection 
of available observations. In the figures, we show the calculated results using 
the default model parameters by the solid black lines, and their estimated 
systematic uncertainties by gray bands. The sensitivity of the results on the 
alternative setups — different Galaxy halo sizes and an alternative source 
distribution — will be discussed below. 

6.1. Anti-proton spectrum and fraction 

The p spectrum and p/p ratio predicted by the model are found to be 
consistent with measurements by PAMELA and other experiments in Figs. H] 
and [5] respectively. The 'SNR-cloud interaction' class contributes in addi- 
tion to the 'SNR-propagation' component to the p intensity with a similar 
spectral shape. Using the latest measurements of p spectrum by PAMELA, 
especially above 10 GeV, it is possible to obtain a crude upper-limit on the 
total number of SNRs interacting with molecular clouds in the Galaxy by 
requiring that the total model p spectrum does not overshoot the PAMELA 
data points and their errorbars. Given the assumed source distribution, un- 
derlying proton spectrum, and typical n x W p (see Section H~2j) . this upper- 
limit is estimated to be ~ 200. In Figs. H] and [51 the maximum possible 
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Figure 4 Measured CR anti-proton spectra compared with the default model. 
Anti-protons from the 'SNR-propagation' class is shown by the solid line. 
The maximum possible contribution of p from the SNR-cloud interaction 
sites (upper-limit) is shown by the dash-dotted line. The dashed line is the 
sum of the two components. 

contributions from the interaction sites are shown for illustration. Of course, 
the typical value of n x W p adopted is derived from a rather small sample 
observed by Fermi, meaning a possible significant systematic error on the 
upper-limit for the total number of these objects in our Galaxy. However, 
their maximum contribution of secondary CR estimated basing on p data 
is independent of the exact value of n x W p used, and should be reasonably 
robust. Future high-precision measurements of the p spectrum and p/p ratio, 
as well as gamma-ray observations of more SNR-cloud interaction systems, 
are expected to provide further constraints on the contributions from these 
secondary CR injectors. 
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Figure 5 Antiproton-to-proton ratio calculated for the default model com- 
pared with measurements. The solid line shows the ratio from the 'SNR- 
propagation' class only. The dash line shows the ratio including also the 
maximum contribution from the SNR-cloud interaction sites (dash-dotted 
line in Fig. @J . 



6.2. Electron + positron (lepton) spectrum 

We see in Fig. [6] that our model agrees with the observed lepton spec- 
trum reasonably well above ~ 20 GeV, within the estimated systematic un- 
certainty. A majority of the CR lepton flux is made up by the SNR-injected 
primary e~ which, however, comes short of reproducing the observation data 
at higher energies; in particular, the points measured by Perm, and 



PAMELA [Hj. The PWN- injected leptons fill up this gap and dominate in 
the higher energies, bringing the total spectrum back to agreement with the 
observed flux. The key for the reproduction of the lepton data is the break 
energy of the PWN injection spectrum at a few hundred GeV. Contribu- 
tion from the SNR-cloud interaction sites is constrained by the upper-limit 
estimated above from the p spectrum, and is found to be almost negligible 
relative to the total spectrum at all energies. The maximum secondary lep- 
ton flux from the SNR-cloud interaction sites is also shown in Fig. [6j In 
terms of the total lepton energy budget in the Galaxy above 1 GeV, contri- 
butions from 'PWN-propagation' class and SNR-cloud interaction sites make 
up for only about 2 %, while the 'SNR-propagation' class is responsible for 
the remaining 98 %. 
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Below ~ 20 GeV, discrepancy is found between the measured lepton 
spectrum and the model. CR spectra in this energy regime are known to be 



strongly affected by the heliospheric modulation [e.g., |87|, and the inaccuracy 



of the force-field approximation applied [88| can be partially responsible for 
the disagreement. We do not observe, however, such a large discrepancy 
with observation in the CR p and p spectra, which would have been affected 
by a similar modulation effect. Although we can not rule out completely 
the modulation effect as the cause of this discrepancy, we will explore an 
alternative possibility later. 

6. 3. Positron fraction 

One general feature predicted by our model is an enhancement of the 
positron fraction above ~ 10 GeV (Fig. [7]), as was recently observed by the 
PAMELA experiment. This enhancement is found to be closely related to the 
hard broken power-law spectrum of the PWN-injected leptons. Without the 
e + contribution from the Galactic PWN ensemble, the 'SNR-propagation' 
source class predicts a strictly decreasing positron fraction with energy in 
contradiction to the PAMELA points. Injection of secondary e + from the 
SNR-cloud interaction sites may boost the positron fraction by < 15 % for all 
energies, but is obviously insufficient to reproduce the 'enhancement' implied 
by the PAMELA data. 

The fractions measured by PAMELA in the 1 — 4 GeV energy range is 
systematically lower than previous measurements. One possible reason is 
a charge sign- dependent heliospheric modulation effect. It has long been 
proposed that the solar modulation of the Galactic CR propagating to the 
Earth depends on the particle charge and mass. The dependency changes 



significantly as the polarity of the solar magnetic field switches [e.g., |89|, [90, 



91| . A significant drop of the positron ratio is indeed predicted during the first 



few years of the PAMELA mission time when the Sun is in its A minimum 



phase [92|. The AMS — 02 team has studied the energy dependence of the 



solar modulation to predict the positron ratio to be reduced significantly 



below ~ 5 GeV |93|, [94J. We reserve a more careful treatment of the solar 
modulation effect for future study. 

On the other hand, the large systematic uncertainty estimated above 
100 GeV mainly arises from the poorly constrained break energy of the 
steady-state PWN injection spectrum, which is seen to vary roughly in the 
broad range of 100 — 500 GeV from gamma-ray observations. 
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Figure 6 Measured CR lepton spectrum compared with the default model 
(thick solid). The thin solid line shows the LIS without solar modulation. 
Contributions from the 'SNR-propagation' and 'PWN-propagation' source 
classes are shown by the dashed and dotted lines respectively. Note that the 
error bars for the PAMELA points are statistical only. The dash-dotted line 
show the maximum contribution from the SNR-cloud interaction sites. The 
gray band shows the estimated systematic uncertainty of the total spectrum 
(see section O for detail). 
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Figure 7 Measured positron fraction compared with that calculated using 
the default model. The fraction calculated for the 'SNR-propagation' class 
alone is shown by the dashed line. The solid line shows the total of the 
'SNR-propagation' and 'PWN-propagation' components, with the gray band 
representing the associated systematic uncertainty. The dash-dotted line 
adds the maximum contribution from the SNR-cloud interaction sites (i.e. 
the dash-dotted line in Fig. [6]) to the solid line, showing the possible effect of 
secondary CR injection from these objects on the observed positron fraction. 
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We found that positrons injected by the Galactic ensemble of PWNe play 
a central role in reproducing the observed positron fraction enhancement. 
There is however a possibility that only PWNe lying within a few 100 pc 
from Earth are significant in contributing to the e + flux. We use the prop- 
agation code described in the Appendix to estimate this local contribution 
from sources with distance-to-Earth smaller than 300 pc, with the same con- 
tinuous source distribution adopted in our default GALPROP model. The 
calculation shows that these nearby PWNe make up less than 10% of the 
total calculated e + spectrum from the 'PWN-propagation' source class in 
the PAMELA energy band. In the steady-state picture of Galactic CR in- 
jection studied in this work, PWNe beyond the local space explain most of 
the observed e + spectrum and fraction. This does not, however, rule out the 
possibility that a small fraction (up to ~ 30 %) of the enhancement observed 
by PAMELA is due to non- steady- state injection from a few nearby PWNe, 



as already suggested in some recent papers [e.g., [16]. 



6.4- A possible low-energy electron component 

A possible cause for the low-energy discrepancy of the lepton spectrum 
is that a class of CR sources (other than the presumed Galactic ensemble of 
SNRs and PWNe) may exist within a few 100 pc of Earth that contribute 
predominately low-energy e~, but not e + . We note that Fermi has found 



gamma-ray emission from X-ray binaries |95|, |96J and Suzaku has found 
gamma-ray pulsar-like periodic hard X-ray emission from a white dwarf 97 
suggesting acceleration of electrons in its magnetosphere. Assuming the ex- 
istence of such low-energy e~ sources, we add an ad-hoc power-law spectrum 
with an exponential cutoff to the calculated lepton spectrum as shown in 
the upper panel of Fig. [HJ The positron fraction is naturally reduced accord- 
ingly as shown in the bottom panel. Since CR e~ with energy ~ 10 GeV 
propagate only a few 100 pc over 1 Myr, these speculative sources should be 
relatively close to Earth. However, it is in order to note that this low-energy 
discrepancy is a result of one of our assumptions that the injected protons 
and primary e~ from SNRs share a common injection spectral shape (for en- 
ergies below E mSLX of e - ), which is supported by the standard DSA theory as 
well as on general physics ground but have no direct observational evidence 
so far. 
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Figure 8 Top panel: Measured CR lepton spectrum compared with the de- 
fault model (thick solid) including a possible new electron component in the 
low energy band (long dash). Other lines are the identical to those in Fig. |6j 
Bottom panel: Measured positron fraction compared with that calculated 
using the default model. The dashed line is the same as the solid line in 
Fig. El while the solid line now shows the modified positron fraction taking 
into account the assumed new electron component. The systematic uncer- 
tainty associated with the modified ratio is displayed by the gray band. The 
dash-dotted line shows the positron fraction when the maximum contribution 
from the SNR-cloud interaction sites is also included. 
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6.5. Gamma-ray spectrum 

We compare in Fig. [9] the gamma-ray spectrum calculated by our model 



for the |6| > 10° region in the Galaxy with that measured by Fermi [98 
which serves as an additional consistency check. The contributions from the 
3 source classes described in section [2] are shown separately. The total spec- 
trum integrati ng oy er all source classes, gamma-ray point sources from the 1st 



Fermi Catalog [100] , and the isotropic background (extragalactic diffuse emis- 
sion and noise induced by CRs) comes out consistent with the data within 
error barsj^l The majority of the observed flux is dominated by the diffuse 
gamma-rays produced by SNR-injected CRs along their propagation path, 
the isotropic background, and the point sources. Diffuse emission produced 
by the PWN-injected leptons are relatively small. Within the upper-limit set 
for the total number of SNRs interacting clouds using high-energy p data, 
gamma-rays emission from the 'SNR-cloud interaction' class is found to be 
insignificant in this region of the sky (the maximum contribution from the 
interaction sites are shown in the figure for illustration). 

6. 6. Effects of modified source distribution model and Galaxy halo heights 

Calculations performed for the default model are repeated for the three 
alternative models - one with an alternative source distribution (see sec- 
tion [3]), and the others with different Galaxy halo sizes. The impacts of the 
alternative setups on our results are summarized below: 

• To re-fit to the CR measurement data, the injection luminosities of 
primary CRs for both the 'SNR-propagation' and 'PWN-propagation' 
classes have to be slightly varied against those adopted for the default 
setup. The required adjustment in the injection luminosities, as well as 
the resulting total CR energy budget in the Galaxy, are summarized in 
Table HJ Most of the deviations are found to stay within ±20 % from 
the default model. 

• For models with different Galaxy halo sizes, the diffusion coefficient 
has to be modified to reproduce the observed B/C ratio. The change 
made is shown in Table [TJ 



6 It is not surprising that the gamma-ray spectrum our model predicts is consistent with 
the Fermi observation because the GALPROP proton injection spectrum has already been 
tuned to reproduce the Galactic gamma-ray spectrum. More rigorous tests will be done 
when high quality analyses of large-scale Galactic diffuse emission become available. 
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Figure 9 Gamma-ray spectrum observed by Fermi from the |6| > 10° region 
(points with error bars) compared with the default model. The total spec- 
trum (solid black) includes the emission produced by SNR-injected CRs (red 
circle) and PWN-injected CRs (blue square) during their propagation, the 
emission from the SNR-cloud interaction sites (upperlimit) through pion- 
decay (green triangle), the isotropic background (cyan diamond), and the 
sum over the 1FGL Fermi catalog point-sources (brown inverted triangle). 
The emission represented by the red circle consists of pion-decay (dashed 
red), IC (dot-dashed red) and bremsstrahlung (dotted red) contributions. 
The emission produced by PWN-injected CRs consists predominantly of IC. 
The Fermi data points are taken from Abdo et al. j98|. 
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Figure 10 Comparison of positron fraction calculated using models with dif- 
ferent assumed Galaxy halo heights. The mean spectrum is used here for the 
PWN-propagation CR injection class for all models, without showing their 
associated systematic uncertainty for clarity reason. The possible contribu- 
tion from the SNR-cloud interaction sites is not included here. 
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• Our results are found to be essentially insensitive to the change of the 
assumed source distribution models, besides the slight changes in CR 
injection luminosities needed. 

• We notice that modification of the assumed Galaxy halo heights have 
subtle influence on the calculated positron fraction (Fig. [TUl) . A smaller 
halo size is found to result in a slightly larger predicted ratio, and vice 
versa. The energy dependence of the ratio is also somewhat affected. 
With the consideration of the spread among the measurement data as 
well as the systematic uncertainty associated to the calculated fraction 
which stems from the assumed PWN injection spectrum, however, the 
observed change in the results are relatively insignificant for any strong 
conclusion to be drawn. The two models with altered halo heights can 
still be considered to be in agreement with measurements. 

• The gamma-ray spectrum shown in Fig. [9] for the > 10° sky region 
does not alter in any significant way among the models. Models with 
different halo sizes changes the overall normalization of the calculated 
spectrum, but by only less than 10% from the default model. More 
detailed comparison with large-scale gamma-ray data in the future, e.g. 
spectra from the Galactic ridge and the inner-Galaxy region, as well 
as small-scale structures, may make it possible to discriminate against 
the models and help put constraints on important quantities such as 
the Galaxy halo height, and the gas and source distributions in our 
Galaxy. 

At this point, we can argue that our results and conclusions obtained for the 
default model are robust against modifications of the source distribution and 
Galaxy halo heights. 

6. 7. Energy budget of Galactic CR 

With gamma-ray observations of CR sources and measurements of charged 
CR spectra analyzed coherently under a common calculation platform using 
GALPROP, it is possible to extract a self-consistent set of quantities re- 
lated to Galactic CRs, including the injected CR energies from SNRs and 
PWNe, interaction rate of nuclear CRs, and gamma-ray emissivity in the 
Galaxy. These quantities depend mildly on the presumed boundary of the 
propagation volume and the assumed source distribution. The results are 
summarized in Table H] and Table Based on the calculated energetics, 
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Table 4. Energetics of Galactic Cosmic Rays in the Steady State 



Particle T p > 0.1 (1.0) GcV Alt. Model l b Alt. Model 2 C Alt. Model 3 d 
[erg] or [erg s" 1 ] A in % A in % A in % 



K.E. of CRs injected by SNRs integrated in the steady-state Galaxy 



+14.8 % 
+9.3 % 
+ 1.8 % 
+ 1.5 % 
-0.5 % 



Injection luminosity of CRs from SNRs 


p 6.4 x 10 40 (4.2 x 10 40 ) +14.6 % 


+3.0 % 


-6.5 % 


e- (pri) 1.2 x 10 39 (8.0 x 10 38 ) +4.1 % 


-23.3 % 


+21.1 % 


K.E. of CRs injected by PWNe integrated in the 


steady-state Galaxy 


r 


e" 6.6 x 10 51 (4.8 x 10 51 ) +9.2 % 






e+ 6.6 x 10 51 (4.8 x 10 51 ) +9.2 % 






Injection luminosity of CRs from PWNe f 


e" 2.3 x 10 37 (2.2 x 10 37 ) +0.7 % 


-19.0 % 


+ 19.1 % 


e+ 2.3 x 10 37 (2.2 x 10 37 ) +0.7 % 


-19.0 % 


+ 19.1 % 



a The results listed here are extracted from the GALPROP outputs for our models. 
b Model using the alternative source distribution (see Section [3]) 
c Model using a smaller Galactic halo height of Z max = 2 kpc. 

d Model using a larger Galactic halo height of Z max = 10 kpc. 

c Kinetic energy of the accumulated CRs within the assumed boundary of the 
Galaxy. Models with different propagation volumes are not compared. 

Kinetic energy and luminosity correspond to the model using the mean spectrum 
for the steady-state injection (see Section [5|). 



p 7.5 x 10 55 (6.1 x 10 55 ) 

e- (pri) 1.9 x 10 54 (7.4 x 10 53 ) 

e- (sec) 1.7 x 10 53 (5.2 x 10 52 ) 

e+ 6.9 x 10 53 (1.5 x 10 53 ) 

p 4.4 x 10 51 (4.4 x 10 51 
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Table 5. 7r°-decay gamma-ray emissivities in different regions of the 

Galaxy a 



Region 


£ 7 > 0.1 GeV 


E 1 > 1.0 GeV 


Alt. Model I c 


Alt. Model II d 


Alt. Model III 6 




[s^H" 1 ] 




A in % 


A in % 


A in % 


Inner-Galaxy a 


2.7 x 10" 25 


3.2 x 10" 26 


-29.7 % 


-29.3 % 


-13.5 % 


Hi-latitude b 


1.6 x 10" 26 


1.9 x 10~ 27 


+10.2 % 


-5.1 % 


+43.3 % 


Entire Galaxy 


1.8 x 10" 26 


2.1 x 10~ 27 


+9.1 % 


+5.1 % 


+33.1 % 



a fi < 5 kpc, \Z\ < 500 pc 
b |Z| > 500 pc 

c Model using the alternative source distribution (see Section [3j 
d Model with Galaxy halo size of Z ma x = 2 kpc. 
c Model with Galaxy halo size of Z max = 10 kpc. 

f The results listed here are extracted from the GALPROP outputs for our models. 



various useful quantitative information can be extracted as we demonstrate 
below for the default set of model parameters: 

• From the total energy and luminosity given for primary protons in 
Table HJ we can roughly estimate the typical lifetime of CR protons 
(T p > 0.1 GeV) inside the Galaxy as ~ 7.5 x 10 55 /2.0 x 10 48 = 3.8 x 
10 7 yr. Similar estimation for CR leptons injected by PWNe is roughly 
9.1 x 10 6 yr, substantially shorter than the protons. This mean lifetime 
is roughly consistent with the synchrotron loss time-scale of electrons 
with an energy of ~ 20 GeV [e.g., |70j], close to the mean energy of 
leptons injected by the PWNe ensemble in our model. The difference 
in the lifetime becomes larger for CRs with T > 1 GeV: ~ 4.6 x 10 7 yr 
for protons versus ~ 6.9 x 10 6 yr for leptons. 

• Approximating the CR proton spectrum in the Galaxy (T p > 1 GeV) 
by a broken PL shape consistent with LIS (E& = 7 GeV, 71 = 2.0, 
72 = 2.75), we obtain the average kinetic energy per CR (T p > 1 GeV) 
to be 3.5 GeV or 5.6 x 10~ 3 erg. 

• For a typical lifetime of ~ 4 x 10 7 yr, protons cross a column density of 
An x 10 25 cm~ 2 where n is the gas density averaged over the propagation 
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volume bounded by R ma x = 30 kpc and Z max = 4 kpc. Since pp 
inelastic cross-section is ~ 30 mb or ~ 30 x 10~ 27 cm 2 , about 1 x n 
of CRs interact with the ISM gas. Since the mean kinetic energy of 



protons is ~ 3.5 GeV, and n multiplicity is ~ 0.5 [61]. Hence a CR 



proton produces, on average, ~lxti pionic gamma-rays in its life. 



For a total number of target H-atom of 1.1 x 10 67 from H I, H2 and 



H II regions given in Ferriere [101] and the pionic gamma-ray emissivity 
averaged over the default GALPROP cylindrical box, the number of 
pionic gamma-rays {E 1 > 1 GeV) coming out of the Galaxy can be 
estimated as (4.5 x 10" 27 ) x (1.1 x 10 67 ) = 5.0 x 10 40 s" 1 . Since the 
total number of CR protons (T p > 1 GeV) is estimated to be 6.1 x 
10 55 /5.6 x 10~ 3 = 1.1 x 10 58 , the average number of pionic gamma-rays 
(£ 7 > 1 GeV) is (5.0 x 10 40 )/(1.1 x 10 58 ) = 5 x 10~ 18 s" 1 per CR 
proton, or (10 15 ) x (5 x 10~ 18 ) = 5 x 10~ 3 integrated over the typical 
lifetime of one CR proton of ~ 10 15 s. Hence, CR protons spend 
most of their life on average propagating in a low-density medium, 
with n ~ 0.005 cm -3 , whereas the typical density of the 'SNR-cloud 
interaction' sites is ~ 100 cm -3 . 

Because of the constraints given by the gamma-ray observations of the 
CR sources and by the CR proton and lepton flux measurements, these es- 
timations are stable within a ~ 20 % level and provide a coherent picture of 
Galactic CR propagation. 

The overall success of the model in reproducing the observed CR ratios 
and spectra suggests that the majority of the CRs observed at Earth are 
steadily injected from the Galactic ensembles of SNRs and PWNe. The 
consistency between the observed and predicted diffuse gamma-ray spectra 
strengthens this interpretation. 

Dependence of the Galactic C R ene rgy budget on CR diffusion models 



has been studied by Strong et al. [103] within the GALPROP framework. 



They report how the diffusion coefficient varies for the 3 halo sizes assumed 
in the present study (2kpc, 4kpc and lOkpc), for the plain diffusion and 
for the diffusive re-acceleration models, under generic constrains that radio, 



Fermi electron [9[ and Fermi diffuse gamma-ray data [58l . l59l . l98l . 199] to be 
reproduced. The CR injection spectra and source distribution have been 
fixed. Our diffusion coefficients and re-acceleration parameters corresponds 
to their Model2 (z04LMS). 
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By comparing the study by Strong et al. [103| with that presented here, 
we note that the CR injection luminosity is much less dependent on the halo 
size. This means Galactic CR modeling is constrained better by CR data 
and gamma-ray measurements near the region where injection is taking place. 
As the Fermi data accumulate further, we will be able to combine the two 
approach and constrain the diffusion process with measurements. 



7. Summary and Conclusions 

The analysis presented here has brought the recent gamma-ray observa- 
tions of CR acceleration sites by Fermi and ACTs and the recent charged CR 
measurements at Earth to one common platform of GALPROP. We were able 
to do this by assuming all CR species in the Galaxy are in a steady state, on 
which GALPROP is built. CRs are assumed to be injected only from SNRs 
and PWNe with the spectra averaged over their respective evolution history. 
The sources are distributed, approximately proportional to, either distribu- 
tion of pulsars or SNRs in the Galacto-centric cylindrical coordinate system. 
Throughout the analysis we have not used the positron fraction neither ex- 
plicitly nor implicitly to constrain the input parameters for GALPROP. The 
analysis has also remained totally blind to the measured p/p ratio. With 
these assumptions, the model reproduces the positron fraction observed by 
PAMELA reasonably well (Fig. EJ), and the p/p ratio by PAMELA (Fig. E]) 
very well under constraint of the observed total e~ + e + spectrum (Fig. EJ) 
0. We have tested the robustness of the model results against alternative 
assumptions on source distribution as well as the Galaxy halo height, and 
confirmed that the conclusions drawn in this study stay intact. 

We conclude that the energy dependence of the excess in positron fraction 
is mostly explained by the contribution from the Galactic ensemble of PWNe 
added upon the secondary e + produced by SNR-injected CR protons along 
their propagation path in the ISM. This does not, however, exclude contri- 
bution coming from local PWNe up to ~ 30 — 40 % in flux. The observed p 
flux is predominantly attributed to interactions of the primary nuclear CRs 
along their propagation paths, while those potentially injected from dense 
clouds interacting with middle-aged SNRs are estimated to be smaller than 
~ 11 % at most in total flux above 100 MeV. The steady-state contribution 



7 We note that our model inherits the prediction of conventional GALPROP models on 
the p/p ratio except for the contribution from the SNR-cloud interaction. 
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from local SNRs and PWNe (within 300 pc from Earth) is estimated to be 
less than ~ 10 % in the e + spectrum, which is not crucial in any way for 
reproducing the observed positron fraction. 

Our result implies that for a typical n x W p (T p > 1 GeV) of 5 x 
10 51 erg/cm 3 , up to ~ 200 SNRs can be interacting with dense clouds and 
injecting secondary CRs (section E]). The maximum contribution of positrons 
at Earth from these interaction sites is found to be relatively insignificant and 
in itself cannot reproduce the observed positron fraction without including 
injection from the PWNe. 

By assuming a SN rate in the Milky Way of one in 30 years 0, the injection 
luminosity of primary CR protons from SNRs listed in Table HJ and SNRs' 
active life of 5 x 10 4 yr, we can estimate that ~ 1700 SNRs are actively 
injecting CRs into the ISM, and each SNR injects ~ 6.1 x 10 49 erg of CRs 
(E > 0.1 GeV) into the Galactic space in its lifetime. Similarly, if we assume 
a PWN birth rate of one in 50 years and that PWNe remain active for 10 5 yr 
on average, ~ 2000 PWNe are injecting CR leptons, with each PWN injecting 
~ 7.2 x 10 46 erg of CR leptons into the Galactic space. CR protons with 
T p > 1 GeV stay ~ 4.6 x 10 7 yr in the Galaxy with only 0.1% interacting 
with gas while CR leptons with T e > 1 GeV lose energy in ~ 6.9 x 10 6 yr. 

A noticeable discrepancy with data is found in the low energy total lepton 
spectrum below ~ 20 GeV. If the corrections for the solar modulation effect 
are properly approximated to the measured spectra, it may imply that an 
unknown additional class of local CR sources which preferentially inject low- 
energy (E < 20 GeV) electrons is required. 

Higher precision gamma-ray spectra expected from Fermi and CTA, as 
well as precision measurements on the positron frac tion above 100 GeV and 



the lepton spectrum in the TeV range by AMS-02 |106j will allow us to ex- 
tract important information on the evolution of PWNe and SNRs. Fermi 
will soon measure anisotropy in the gamma-ray emissivity at nearby molec- 
ular clouds, and either detect or exclude the proposed GRB scenario for the 



PAMELA e + excess [19|, |56j. Future high precision measurement of p by 
AMS-02 will further constrain contribution from the SNR-cloud interaction 
sites. 



8 Diehl et al. |102| estimate the rate of core-collapse supernovae to be ~ 1.9 ± 1.1 per 
century. The Type la su pernovae rate is less constrained observationally, ranging from 
0.3 — 1.1 per century [e.g. 104 , 105) . These two rates add up to ~ 3 per century. 
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Appendix A. A Supplementary Propagation Scheme 

Models described in this study are primarily built using GALPROP as 
far as possible. Other calculations that cannot be readily implemented in a 
typical GALPROP model are accomplished via a simple propagation code 
separately developed by us, which include the following: (1) the local con- 
tribution of CR positrons from PWNe lying within a given distance from 
Earth; (2) contribution of secondary particles injected from SNRs interact- 
ing with clouds; and (3) the high-energy rollover inherit to the primary e~ 
spectrum injected from SNRs due to spectral cutoff at a time-evolving maxi- 
mum energy. To obtain the CR spectra propagated from a continuous source 
distribution characterized by distance-to-Earth r, source age t age , and the 
active lifetime of each source tnf e (time interval during which CRs are in- 
jected into the ISM from a source), we use the analytic solution of the simple 
diffusion-loss transport equation adapted from Atoyan et al. [42] 0: 

QmjjEj, r, t')b(Ej) -(\r\/r D (E,t>)f dt > (A ^ 



9 This propagation calculation does not include re-acceleration effect and interaction 
involving heavy nuclei which can possibly affect the B/C ratio. 



N(E,r)= [ 
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r D (E,t) «2 



1 - (1 - b Ety~ s 
(1 - 5)6 £t 



(A.2) 



Here b(E) is the energy-loss rate of the CRs, which can be expressed as 
b(E) ~ boE 2 for synchrotron and IC losses of leptons during propagation, 
while loss is considered negligible for nuclei (such that equation IA.1I and IA.2I 
reduce to the usual solution of a simple 1-D diffusion equation and diffusion 
length respectively). E t is the initial injection energy of CRs whose observed 
energy at Earth is E. The effective diffusion length rn(E,t) for a particle 
with observed energy E and propagation time t is given by equation IA.2[ in 
which energy loss effect is accounted for. D(E) oc f3E 5 is the adopted spatial 
diffusion coefficient, where the index S = 1/3 is used in consistence with the 
GALPROP calculations. The CR injection term is given by Q- m j(E, r,t) = 
Q(E,t)n BOC (r,t) where Q(E,t) is the injection spectrum and n soc (r,t) is 
the source distribution in space-time. Spatial distribution of sources is the 
same as the distribution used in the default GALPROP model, while t age of 
the sources are assumed to be distributed uniformly in time starting from 
a maximum of 10 8 yr up to 5 x 10 4 yr from now. For every source, tuf e is 
assumed to be 5 x 10 4 yr. 

For (1), the injection spectra and luminosities of e + and e~ are di- 
rectly taken from the GALPROP inputs and outputs adopted for the 'PWN- 
propagation' class of the default model (see Tables Q] and HJ). For (2), these are 
calculated for the secondary e + , e~ and p according to the description in Sec- 
tion 14.21 Tertiary CRs produced through hadronic interactions during prop- 
agation of p in the ISM are neglected. For (3), a time- dependent exponen- 
tial cut-off is applied to the injection term for the SNR-injected primary e~ 
(Q(E,t) oc e~ E ^ Eraa - x ^) to incorporate the time-evolving maximum electron 
energy E max (t) given by equation [2j This results into a smooth high-energy 
turnover of the propagated primary e~ spectrum at Earth, which commits at 
around 1 TeV. This calculated turnover is then applied to the GALPROP- 
calculated primary e~ spectrum from the 'SNR-propagation' source class. 
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